arXiv:1507.02434vl [astro-ph.SR] 9 Jul 2015 


Spectral Variations of Of?p Oblique Magnetic Rotator Candidates 

in the Magellanic Clouds 

Nolan R. Walborn 

Space Telescope Science Institute^, 3700 San Martin Drive, Baltimore, MD 21218 

walbornOstsci.edu 
Nidia I. Morrell 

Las Campanas Observatory, Carnegie Observatories, Casilla 601, La Serena, Chile 

nmorrellOlco.cl 
Yael Naze^ 

GAPHE, Departement AGO, Universite de Liege, Allee du 6 Aout 19c, Bat. BSC, 

B4000-Liege, Belgium 

nazeQastro.ulg.ac.be 

Gregg A. Wade 

Department of Physics, Royal Military College of Canada, PO Box 17000 Station Forces, 

Kingston, ON, Canada K7K 7B4 

wade-g@rmc.ca 

Stefano Bagnulo 

Armagh Observatory, College Hill, Armagh BT61 9DG, UK 

sbaSarm.ac.uk 
Rodolfo H. Barba^ 

Departamento de Fisica y Astronomia, Universidad de La Serena, Cisternas 1200 Norte, 

La Serena, Chile 

rbarbaOdfuls.cl 

Jesiis Maiz Apellaniz^ 

Centro de Astrobiologia, CSIC-INTA, Campus ESAC, Apartado Postal 78, E-28691 

Villanueva de la Canada, Madrid, Spain 


- 2 - 


jmaizOcab.inta-csic.es 
Ian D. Howarth"^ 

Department of Physics and Astronomy, University College London, Gower Street, London 

WCIE 6BT, UK 

idhOstar.ucl.ac.uk 

Christopher J. Evans^ 

UK Astronomy Technology Centre, Royal Observatory Edinburgh, Blackford Hill, 

Edinburgh EH9 3H.J, UK 

Chris .evansQstfc.ac.uk 

Alfredo Sota 

Instituto de Astrofisica de Andalucia-CSIC, Glorieta de la Astronomia s/n, 18008 

Granada, Spain 

sota@iaa.es 

ABSTRACT 

Optical spectroscopic monitoring has been condncted of two O stars in the 
Small and one in the Large Magellanic Clond, the spectral characteristics of which 
place them in the Of?p category, which has been established in the Galaxy to 
consist of obliqne magnetic rotators. All of these Magellanic stars show sys¬ 
tematic spectral variations typical of the Of?p class, farther strengthening their 
magnetic candidacy to the point of virtnal certainty. The spectral variations 
are related to photometric variations derived from OGLE data by Naze et ah 
(2015) in a parallel stndy, which yields rotational periods for two of them. Now 
circular spectropolarimetry is required to measure their helds, and ultraviolet 
spectroscopy to further characterize their low-metallicity, magnetically conhned 
winds, in support of hydrodynamical analyses. 
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Subject headings: Magellanic Clouds — stars: early-type — stars: emission-line 
— stars: magnetic fields — stars: rotation — stars: variables 


1. Introduction 

The Of?p class in the Galaxy currently consists of five objects (Walborn et ah 2010). 
The sometimes maligned question mark in their classihcation, introduced by Walborn (1972, 
1973), was intended to denote doubt that they are normal Of supergiants, as the latter 
category was being interpreted at the time, because of several relatively subtle spectral 
peculiarities. One of them is emission-line prohles narrower than those of the absorption lines, 
suggesting a circumstellar origin of the former in the peculiar stars, rather than photospheric 
as in the normal ones. That morphological distinction was amply vindicated decades later, 
when 1-20 kG magnetic fields were discovered in all of the Of?p stars (Donati et al. 2006; 
Martins et al. 2010; Wade et al. 2011, 2012a, 2012b, 2015; Hubrig et al. 2008, 2011, 
2013, 2015). Because of obliquity between their magnetic and rotational axes, they display 
spectacular variations during their rotational cycles, which range from about a week to more 
than 50 years among the Galactic objects, likely determined by spindown as a function of 
field strength and age, due to angular momentum loss produced by coupling of the magnetic 
helds to the winds (“magnetic braking”). That is, they are higher-mass analogues of the Ap 
and Bp stars, with even more spectacular behaviors because the helds channel the strong 
stellar winds from opposite hemispheres into a collision at the magnetic equator, forming a 
thin disk that presents varying aspects. While on the one hand there is a phenomenological 
unity among the class because of the basic mechanism, on the other hand the different 
combinations of multiple parameters also lead to a remarkable diversity of detail. 

This diversity among the small number of Galactic Of?p stars is not merely amusing; 
the ranges among their individual parameters provide “separated-variable experiments” that 
illuminate the entire class. Glearly further examples of the phenomenon are of considerable 
interest from that perspective. The hrst Of?p star in the Magellanic Glouds was recognized 
in the SMG by Walborn et al. (2000-still before the discovery of their physical nature!); 
two more in each Gloud have been found since. They offer the interesting opportunity to 
investigate magnetic helds and winds of massive stars at lower metallicities. Here we present 
spectroscopic monitoring of the three objects known for long enough to accomplish that; 
two more have been discovered very recently by Massey et al. (2014) and are not further 
discussed in this paper, although they are included in the companion photometric paper by 
Naze et al. (2015). 
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2. Observations 

Spectroscopic observations of AV220 and 2dFS936 in the SMC, and BI57 in the LMC 
(designations explained below), have been obtained with the Boiler & Chivens spectrograph 
attached to the 2.5ni du Pont telescope at Las Campanas Observatory. The Marconi No. 1 
2048x515 13.5yu-pixel detector was in use with that instrument. We used a 1200 lines mm“^ 
grating centered at 4700 A and had the slit width set to 150 microns, corresponding to 
1726 on the sky and 1.75 pixels on the detector. This instrumental conhguration produces a 
resolution of ~1.7 A as measured from the FWHM of the comparison lines. The typical peak 
signal-to-noise ratio (S/N) per 2-pixel resolution element ranges from 150 to 200, with just 
a few cases signihcantly below or above. All of these data were obtained under the technical 
auspices of the Galactic 0-Star Spectroscopic Survey (GOSSS; Maiz Apellaniz et al. 2011), 
and they have been incorporated into the GOSSS database. 

Dome flats and bias frames were obtained during the afternoon prior to each observing 
night as well as twilight flats at sunset, and He-Ne-Ar comparison-lamp exposures were 
recorded before or after each target was observed. Reductions were carried out using standard 
IRAF^ longslit routines. Flux calibration was derived from a set of standard stars observed 
during the same nights, and the data were then normalized for presentation here. 

One observation of AV220, one of 2dFS936, and two of BI57 were obtained with the 
Magellan Echellette (MagE; Marshall et al. 2008) at the Glay 6.5m telescope, with a 1" x 5" 
slit, yielding a resolution ranging from 0.7 A at 3200 A to 1.7 A at 10,000 A, as measured 
from the FWHM of the Th-Ar comparison lines. The peak S/N per 3-pixel resolution 
element is 95 for AV 220, 215 for 2dFS 936, and 140 for BI 57. Reductions were done with 
a combination of the IRAF mtools package (originally designed by Jack Baldwin to process 
Magellan Echelle-MIKE-data) and standard IRAF echelle routines, as described by Massey 
et al. (2012). 

We also present here one observation of 2dFS 936 and one of BI 57 obtained by IDH and 
GJE with the 2dF instrument attached to the Anglo-Australian Telescope (see Lewis et al. 
2002; Evans et al. 2004). One spectrogram of 2dFS 936 was obtained by IDH with the ANU 
2.3m telescope and Dual-Beam Spectrograph (DBS) at Siding Spring Observatory, described 
by Rodgers et al. (1988), using 1200 lines mm“^ gratings and E2V detectors windowed to 
2048x512 13.5/i pixels, which gave a nominal resolution of ~1 A. Also discussed here for 
AV220 are one observation from GASPEG at the 3.6m on ESO/La Silla (Walborn et al. 


^IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associa¬ 
tion of Universities for Research in Astronomy (AURA), Inc., under cooperative agreement with the National 
Science Foundation. 
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2000), and data obtained with the Sm VLT/FLAMES on ESO/Paranal by CJE (LR02/LR03 
conhgurations as described by Evans et ah 2011). 


3. Results 

All of the Galactic Of?p stars display distinctive, periodic spectral variations to various 
degrees, which are reproducible across the rotational cycles for as long as they have been 
observed (Naze et ah 2001, 2008, 2010; Walborn et ah 2004; Howarth et ah 2007). The 
phases of their magnetic extrema correlate with those of the spectral variations, which may be 
described as maxima and minima in terms of the emission-line intensities. These variations 
include the apparent spectral types, with the earlier types due to emission hlling at the 
He I lines, so that the minima are likely more representative of the actual types. In some 
objects, the C III A4650 emission lines, whose comparability to the N III A4640 was one of 
the original criteria for the class, disappear abruptly and completely, as if in an occultation. 
In some objects, Ha changes from a strong emission line to an absorption, and comparable 
prohle changes occur in He H A4686. Higher Balmer series members develop variable, narrow 
P Cygni prohles within the broad stellar absorption lines. Several Of?p stars display very 
small-amplitude light variations that also correlate with the projected magnetic held and 
emission lines, because scattering electrons believed to produce the photometric maxima are 
concentrated over the magnetic poles, corresponding to face-on aspects of the disks which 
maximize the emission-line strengths (Walborn et al. 2004; Howarth et al. 2007; Wade et 
al. 2011, 2012a, 2015). 

Here we report precisely these kinds of variations in the spectra of all three of the 
MC objects discussed, as well as correlations with their light variations, thus essentially 
establishing their Oblique Magnetic Rotator (OMR) nature in advance of any magnetic 
observations. Once again, the ongoing power of morphology to discover and characterize 
phenomena prior to physical measurement and analysis is demonstrated. 


3.1. AV 220 in the SMC 

AV 220 is located just north of NGC 346, the largest H H region in the SMG, but it is a 
held star sufficiently removed to be unaffected by nebular emission lines. A chart is provided 
by Azzopardi et al. (1975, Ghart 7). AV 220 has V = 14.50 and My = —4.9 (for R = 3 and 
a distance modulus of 19.1), the latter far fainter than for a supergiant (Walborn 1973). In 
fact, since the distances of most Galactic Of?p stars are uncertain, these MG objects provide 



Table 1. 


Summary of spectroscopic observations of AV 220 ordered by epoch 


Date 

HJD-2450000.0 

SpT 

EW A4686 

EW H/5 

Instrnment 

20 Nov 1997 

0772.5 

06.5 

-1.53±0.05 

-2.50± 0.07 

CASPEC 

27 Sep/04 Oct 2004“ 

3278.5 

06.5 

-2.00T0.10 

-2.75±0.10 

FLAMES 

03 Dec 2011 

5898.6 

07 

-1.15±0.05 

-1.20T0.10 

B&C 

06 May 2012 

6053.9 

07.5 

-0.86±0.10 

-0.96±0.10 

B&C 

10 May 2012 

6057.9 

07.5 

-0.85±0.10 

-0.75±0.10 

B&C 

16 Mar 2013 

6368.5 

06.5 

-1.15±0.05 

-1.70±0.10 

B&C 

29 May 2013 

6442.9 

06 

-1.35±0.05 

-1.80±0.10 

B&C 

16 Jnn 2014 

6824.9 

05.5 

-2.02±0.05 

-2.57±0.10 

B&C 

18 Jnn 2014 

6826.9 

05.5 

-2.00±0.05 

-2.73±0.10 

B&C 

28 Jnl 2014 

6866.8 

05.5 

-2.10T0.05 

-2.89±0.10 

B&C 

09 Jan 2015 

7031.5 

06 

-2.00±0.05 

-2.90±0.10 

MagE 


'^This observation combines 4 individual spectrograms obtained 2004 Sept. 27 & 28 
(covering AA3980-4560) and Oct. 1 & 4 (AA4510-5050); the Julian Dates have been 
averaged. 








Wavelength (A) 

Fig. 1.— Spectroscopic monitoring of AV 220 in the SMC. Since there is no period for this 
object as yet, the data are plotted in chronological order. Nevertheless, large variations 
typical of Of?p stars in spectral type (He I/He H absorption-line ratios) and emission-line 
intensities are seen (cf. Table 1). Wavelengths of the identihed lines may be found in the 
spectral atlas of Sota et ah (2011). 




I magnitude (OGLE III —IV) 



2000 3000 4000 5000 6000 7000 

HJD-2450000. 


Fig. 2.— Secular plot of the He II A4686 and H/5 emission-line EWs in AV 220, compared 
with its lightcurve from Naze et al. (2015), where details of the photometry can be found. 
See text for discussion. 
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valuable confirmation of their spectroscopic nonsupergiant interpretation. 

Table 1 specifies and Figure 1 displays the 11 available spectrograms of AV 220. The 
1997 CASPEC data (Walborn et al. 2000, Fig. 6) comprise the hrst Of?p spectrum recog¬ 
nized in the Magellanic Clouds. (Note also in Fig. 5 of that reference the high-quality HST 
UV observation made 24 May 1999, which displays very peculiar stellar-wind prohles and 
constitutes the only UV data for an MC Of?p star to date.) 

The apparent spectral type of AV 220 ranges from 05.5 to 07.5, i.e., the He I A4471 /He H A4541 
ratio strongly inverts about the unit value at 07. The Balmer lines are highly variable, as 
best seen at Hy in the hgure, with central or redshifted emission ranging from very weak 
to very strong in correlation with the spectral type (which, again, varies due to emission in 
He I A4471; the weaker He I lines are too strong for the earlier types, which is another char¬ 
acteristic of Of?p spectra no doubt related to the illusory maximum types). The equivalent 
widths (EWs) of the He H A4686 and H/3 emission lines (Table 1) vary correlatively with 
the spectral types. On the other hand, any variations in the C HI/N HI emission ratio are 
at the limit of significance in these data; certainly the C HI never disappears. The essential 
constancy of the spectrum from mid-June through July 2014 is indicative of a long period. 

Unfortunately, despite an extensive database, clear variability, and exhaustive searches, 
the photometry has thus far not yielded a dehnitive period for this star (Naze et al. 2015). 
Therefore, it is not yet possible to relate the substantial spectroscopic variations to its 
(presumed) rotational phases. Nevertheless, a secular plot of its lightcurve and emission-line 
EWs is presented in Figure 2, which shows an apparent correlation in the expected sense (i.e., 
weaker lines at a light minimum and viceversa) during the last ~1000 days of overlapping 
coverage. A timescale of roughly 1500 days could describe this segment. However, the earlier 
photometry is not consistent with that, nor is the (multiple) spectroscopic observation at 
(average) HJD 2453278.5 (Table 1). Evidently, another phenomenon in addition to rotation 
complicates the variations of AV 220. 


3.2. 2dFS 936 in the SMC 

This object was first reported as Of?p by Massey & Duffy (2001, Fig. 8), who provide 
a chart (their Fig. 6); it was subsequently observed and discussed by Evans et al. (2004, 
Fig. 17), whence the designation adopted here. 2dFS 936 has V = 14.02 and My = —5.6, 
again fainter than a supergiant. 

Our observational and spectroscopic data for 2dFS 936, consisting of 9 spectrograms, 
are presented in Table 2 and Figure 3; they are ordered by phase as discussed below. The 



Normalized Flux 



Wavelength (A) 

Fig. 3.— Spectroscopic monitoring of 2dFS 936 in the SMC. This star has a well determined 
photometric period of 1370±30 days so the spectroscopic data are ordered by phase (cf. 
Table 2 and discussion in the text). 




I magnitude (OGLE II-III-IV) 
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2dFS 936 



0 0.5 1 
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Fig. 4.— Phase plot of the He II A4686 and H/3 emission-line EWs in 2dFS 936, compared 
with its lightcnrve from Naze et ah (2015). 
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apparent spectral type ranges from 06.5 to 07.5, with the He I/He H classihcation ratio 
again inverting through unity. The Balmer emission at both H/3 and Hy varies signihcantly 
in intensity correlatively with the spectral types, but not in position, appearing always 
redshifted. The C HI/N HI emission-line ratio is much smaller than in AV 220 and is 
essentially constant. In view of the very long photometric period discussed next, the entire 
2011-2014 spectroscopic sequence corresponds to a single rotational cycle. 

The EWs of the He H A4686 and H/3 emission lines are listed in Table 2 and plotted 
in Figure 4 in comparison with the photometric lightcurve from Naze et ah (2015), which 
corresponds to a well determined period of 1370±30 days. Although sparse, the spectro¬ 
scopic data show a tendency to track the lightcurve in the sense expected from Galactic 
counterparts, i.e., stronger emission lines near the light maximum. The large ratio of the 
He H to the Balmer line is noteworthy. 


3.3. BI 57 in the LMC 

The Of?p nature of this object was detected in an extensive, unpublished LMC survey 
by IDH. It was subsequently identihed with BI 57 (Brunet, Imbert et al. 1975), who provide 
a chart in their Plate 1 (second Hodge-Wright Atlas 28B cutout). It has V = 13.97 and 
My = —5.0 (for a DM of 18.6), similar to the other two stars above and again much fainter 
than a supergiant, as expected and now confirmed for Of?p stars. 

The data for 8 spectroscopic observations of this star are presented in Table 3 and 
Figure 5; the one at phase 0.10 of the 400d period is from the discovery survey. The short 
spectral-type range is always later than 07, i.e., 07.5 to 08, which may seem small but 
is morphologically (and presumably physically) signihcant. In this spectrum the variable 
Balmer emission is always either central or blueshifted, while there is a small but dehnite 
variation in the C HI/N HI emission ratio (compare A4650 to A4634). On the other hand, the 
relatively weak He H A4686 emission is essentially constant whereas H/3 is stronger and more 
variable (cf. EWs in Table 3), the opposite of the case for 2dFS 936. In some spectra there 
is a weak absorption feature longward of the He H A4686 emission, creating the appearance 
of an inverse P Cygni profile. 

Apart from the absolute magnitude and Of?p variations, there is no question of this 
spectrum arising from a normal (or even Ofc-Walborn et al. 2010) LMC supergiant. At 
this late O type, there would be very strong Si IV absorption lines flanking Hh that are 
not present, and the Balmer absorption itself is too strong for a supergiant. Quantitative 
analyses of the high-resolution data for all three objects discussed here will be performed in 



Table 2. Summary of spectroscopic observations of 2dFS 936 ordered by phase 


Date 

HJD-2450000.0 

Phase“ 

SpT 

EW A4686 

EW H/3 

Instrument 

15 Jun 2014 

6823.9 

0.07 

07.5 

-2.43±0.10 

-0.65±0.2 

B&C 

17 Jun 2014 

6825.9 

0.07 

07.5 

-2.44±0.10 

-0.62±0.2 

B&C 

28 Jul 2014 

6866.8 

0.10 

07.5 

-2.75±0.10 

-0.78±0.2 

B&C 

03 Sep 2014 

6903.7 

0.12 

07 

-3.06±0.05 

-1.37±0.2 

MagE 

30 Sep 1999 

1452.0 

0.15 

06.5 

-3.40T0.10 

— 

2dF 

03 Dec 2011 

5898.5 

0.39 

06.5 

-3.20T0.10 

-0.99T0.2 

B&C 

05 May 2012 

6052.9 

0.50 

06.5 

-3.31T0.10 

-1.42T0.2 

B&C 

16 Mar 2013 

6367.5 

0.73 

07.5 

-3.00T0.10 

-0.92T0.2 

B&C 

31 Dec 2005 

3736.0 

0.80 

07.5 

-2.60T0.05 

-0.98±0.2 

DBS 


^Based upon P = 1370±30 days and To=2453993.3±1.7 (Table 2 of Naze et al. 2015) 


Table 3. Summary of spectroscopic observations of BI57 ordered by 400d phase 


Date 

HJD-2450000.0 

Phase“ 

Phased 

SpT 

EW A4686 

EW H/3 

Instrument 

03 Dec 2004 

3343.0 

0.10 

0.51 

08 

-0.43±0.05 

-0.98±0.05 

2dF 

09 Jan 2015 

7031.7 

0.32 

0.20 

07.5 

-0.55±0.05 

-2.10±0.05 

MagE 

17 Feb 2013 

6340.5 

0.60 

0.32 

07.5 

-0.50±0.10 

-1.74±0.10 

B&C 

16 Mar 2013 

6367.6 

0.66 

0.35 

07.5-08 

-0.51±0.05 

-1.56±0.05 

B&C 

13 Feb 2012 

5970.6 

0.67 

0.85 

08 

-0.50±0.05 

-1.10±0.05 

MagE 

25 Apr 2014 

6773.5 

0.68 

0.87 

08 

-0.52T0.05 

-0.97±0.05 

B&C 

04 May 2012 

6052.5 

0.88 

0.95 

08 

-0.53±0.05 

-1.26±0.05 

B&C 

28 Jul2014 

6866.9 

0.91 

0.99 

08 

-0.60±0.05 

-1.42±0.05 

B&C 


^Based upon Pi = 400±3.5 days and To=2455302.2±1.2 (Table 2 of Naze et al. 2015) 
'^Based upon P 2 = 787±14 days and the same Tq (Naze et ah 2015) 
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Fig. 5.— Spectroscopic monitoring of BI 57 in the LMC. The observations are ordered by 
phase corresponding to a period of 400±3.5 days (cf. Table 3 and discussion in the text). 
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due course. 

Although not well distributed, the phase behavior of the He II A4686 (constant) and 
H/5 emission EWs is compared with both the 400±3.5 and 787±14-day lightcurves found by 
Naze et al. (2015) in Figure 6. While there is some semblance of a correlation with H/9 on 
400d, it is not entirely satisfying. The photometric searches found a second value of 787d 
with comparable probability from different methods and epoch ranges, i.e., essentially twice 
the first value. Here two light maxima appear near phases 0.2 and 0.8, unequal particularly 
in their widths. The correlation of the H/3 EWs with the first maximum is remarkable, far 
better than on the 400d period; unfortunately, the second light maximum has inadequate 
spectroscopic coverage so any correlation is unclear. If the 787d period is correct, then 
the 400d period is of course a harmonic of it; and it is likely that both magnetic poles 
are visible during the rotational cycle (as in the Galactic magnetic O stars HD 57682 and 
CPD —28° 2561; Grunhut et al. 2012, Wade et al. 2015, respectively), albeit unequally in 
this case. 


4. Discussion and Outlook 
4.1. Magnetic Observations 

From detailed spectroscopic monitoring of three MG Of?p objects, we have found dis¬ 
tinctive variations analogous to those displayed by their Galactic counterparts, all of which 
have been established as OMRs. Thus, with high conhdence these MG objects also belong 
to that category and they are identihed as the hrst extragalactic magnetic O stars in ad¬ 
vance of any magnetic observations. Nevertheless, it is of course essential to measure the 
strengths, orientations, and (aspect) variations of their helds directly, in order to support 
physical modeling; in particular, the inclination of the rotational axis and obliquity of the 
magnetic axis must be known to model the spectroscopic and light variations. The only 
current instrument capable of such measurements is FORS at the ESO VLT (Bagnulo et al. 
in preparation). Since we have now derived rotational ephemerides for two of them (Naze 
et al. 2015), we can predict the most favorable epochs for those observations, i.e., at the 
photometric and spectroscopic maxima. 


4.2. Ultraviolet Observations 

The UV stellar-wind prohles and variations, especially in their blueshifted absorption 
troughs, provide vital information about the structure of the OMR magnetospheres. As 



16 


discussed above, the field channels the winds from opposite hemispheres into a collision at 
the magnetic equator, where high-energy phenomena are produced. As the magnetosphere 
rotates rigidly and obliquely, the variable wind profiles diagnose the structure and flows in 
different directions. Indeed, three Galactic Of?p stars observed at multiple phases with HST 
display very peculiar and variable UV wind profiles, which can be understood in terms of 
magnetohydrodynamic models (Marcolino et ah 2012, 2013; Naze et ah submitted). A 
further Galactic object is currently under observation (Petit et al. in preparation). As also 
mentioned above, a single HST observation of AV 220 (Walborn et al. 2000) likewise exhibits 
peculiar wind profiles, but it is the only extant UV observation of an MG Of?p star. Now 
it is essential to obtain UV observations at the magnetic and optical extrema of the two 
objects with periods discussed here, and to pursue the period of the third for that purpose, 
to enable full characterization of their magnetospheric structures. 
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Fig. 6.— Phase plots for 400d (left) and 787d (right) periods of the He II A4686 and 
emission-line EWs in BI 57, compared with the corresponding lightcurves (Naze et ah (2015). 














